The Poisson ratio of crystalline surfaces 
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Abstract 

A remarkable theoretical prediction for a crystalline (polymerized) surface is 
that its Poisson ratio a is negative. Using a large scale Monte Carlo simulation 
of a simple model of such surfaces we show that this is indeed true. The 
precise numerical value we find is a ~ —0.32 on a 128 2 lattice at bending 
rigidity k = 1.1. This is in excellent agreement with the prediction a = —1/3 
following from the self-consistent screening approximation of Le Doussal and 
Radzihovsky. 

Crystalline surfaces have been studied extensively in recent years. Unlike one- 
dimensional polymers, which are always crumpled, non self-avoiding crystalline sur- 
faces undergo a continuous phase transition from a high temperature crumpled phase 
to a low temperature flat phase [QJ3|. The flat phase is characterized by long-range 
orientational order in the normals to the surface. 

There are several experimental realizations of crystalline surfaces. Some, like sus- 
pended layers of graphite oxide in aqueous solution j^lHj or polymerized adsorbed 
mono-layers, are chemical systems one can synthetize in the laboratory. There are 
also beautiful biological examples of such surfaces 0: the cytoskeleton of erythro- 
cytes (red blood cells) is composed of a network of nodes (actin oligomers) and links 



1 



Table I: The scaling exponents. 
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(spectrin tetramers). A typical skeleton is a triangulated network of roughly 70,000 
plaquettes. 

The crumpling transition and the flat phase of crystalline surfaces have been 
investigated numerically by several authors. The interested reader may consult the 
excellent reviews [TJIH]- Bowick et. al have recently performed a large scale simula- 
tion of a triangulated crystalline surface with bending rigidity using free boundary 
conditions 9 . The largest surface simulated has 32,258 triangles. The equilibrium 
distribution is sampled using a Monte Carlo algorithm with a local Metropolis up- 
date. The action used has a simple Gaussian potential and a bending energy term: 

/3H = ^2(n-fj) 2 - n a -np, (1) 

(ij) (a/3) 

where r*j is the position of node i, n a is the unit normal to triangle a and k is the 
bending rigidity. For k > k c ~ 0.79 the system is in a flat phase and it behaves like 
a membrane with anomalous elasticity. 

The flat phase is characterized by the scaling of the renormalized effective elastic 
constants Xr ~ fj,R ~ q Vu and of the bending rigidity kr ~ q~ v . The exponent rj 
can be determined from the scaling of the height-height correlation function (rj is 
related to the roughness exponent £ = 1 — rj/2). The exponent rj u determines the 
finite size scaling of the mean square phonon fluctuations. 

We summarize the results of |3j in Table [fl These results are compared to the 
analytical predictions obtained from an e = 4 — D expansion (AL) |10Ullj . a large- 
d expansion |11H13| and a self-consistent screening approximation (LR) In 
this notation D represents the dimensionality of the elastic manifold and d is the 
dimensionality of the embedding space (for physical membranes D = 2 and d = 3.) 
Scattering experiments on the red blood cell skeleton give £ ~ 0.65(10) [S]. 

In addition to the anomalous scaling of the coupling constants, one of the most 
dramatic effects of fluctuations on crystalline surfaces is the prediction of a negative 
Poisson ratio a. The Poisson ratio measures the in-plane transverse response of the 
surface when stress is applied in the longitudinal direction. It is defined to be positive 
for matter which shrinks in the ^-direction when stretched in the y-direction. 

Analytical calculations in the context of a self-consistent screening approxima- 
tion (LR) ^1] and an e expansion (AL) ^U] predict that for crystalline surfaces a 
is —1/3 (LR) and —1/5 (AL) respectively. The unusual sign of the Poisson ratio 
is a result of entropic suppression of the height fluctuations in a membrane under 
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stress This effect is clearly demonstrated by crumpling a sheet of paper and 
pulling on two opposite corners: the sheet expands in the direction transverse to 
the applied strain [TH| [T7] . 

In this letter we demonstrate numerically that a crystalline surface, defined by 
Eq. (fljl. indeed has a negative Poisson ratio, and that it agrees with the prediction 
from [Tl]. We compare our result with previous numerical determinations of a |18j . 
We stress that this work extends the analysis of Monte Carlo data collected in 
and we refer to it for the details of the simulations. 

Consider an asymptotically flat elastic surface at thermal equilibrium. In the 
Monge gauge, its behavior is described by the partition function 



Z[<7«] 
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where ? is the intrinsic coordinate, fiQ and Ao are the bare Lame coefficients and kq 
is the bending rigidity. The stress tensor Oij represents an external source linearly 
coupled to the system. The strain tensor Uy, to linear order in u, is related to f by 



s + u + zh 



1 



(diitj + djUi + dihdjh). 



(3) 
(4) 



Here s is the rest (equilibrium) position of the surface, assumed to lie in the x-y 
plane, and h is the height of the surface above the reference plane. 

The Poisson ratio can be defined in terms of correlation functions at zero external 
stress using linear response theory. Considering a diagonal (hydrostatic) stress, the 
derivation is straightforward and leads to 



a 



(u X x 



yy/c 



(5) 



where the subscript c indicates the connected part. Since the correlation functions 
are measured in the limit of zero external stress (cr^ = 0) the system is isotropic in 
the x-y plane and 

(«L)c = { u ly)c- (6) 

For computational purposes it is convenient to express the strain tensor in terms 
of the tangent vectors ti = dir. The index i refers to the intrinsic coordinate system. 
The strain tensor is 



Uij 



dif- djr 



(dir) ■ (djf) = gij - 5, 



(7) 



where gij = dif ■ djf is the induced metric, and we have rescaled the intrinsic 
coordinates so that (gij) = Sij. Substituting Eq. in Eq. © we get 



a 



(d xx 9yy)c 
(9 2 yy)c 



(8) 
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Table II: The Poisson Ratio, 
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In terms of a discretized surface the induced metric g assumes a simple form. Take 
for example the xx component at point ? 



which is the squared length of the "link" between the point at ? and its neighbor in 
the x direction. 

The derivation of a implicitly assumed that the intrinsic coordinate system is 
orthogonal. As our discretized surface is a triangular lattice, we need to transform 
the coordinate system. In fact, we can directly access only the <9jr f , % = 1, 2, 3, 
of Eq. ®, in the three natural directions of the triangular lattice, while Eq. (jHJ) is 
expressed in terms of x and y. If ex, &2, £3 are the basis vectors of a triangular 
lattice, we define the x direction to overlap with ex, and we use y = (£2 + es)/V3. 

At this point we can use Eq. (JBj), isotropy, as a consistency check on our definition 
of the strain fluctuation. In fact we verified that (g xx ) c , (<?ii)c (<?22)c> (533)0 and 
(dyy)c are equal within numerical accuracy. 

One can find an equivalent definition of a in terms of different correlation func- 
tions. Since {g xx g xx ) c - {gxx9 yy )c = 1{g xy g xy ) c HH|, it is easy to show that 



We have verified that this definition gives results consistent with Eq. (JHJ). But it 
must be noted that it is more difficult to use numerical methods to determine a 
from Eq. (JTHJ) as idxy) is much smaller than {g xx )- 

We report in Table Hl| the measured values of the Poisson ratio for various sizes of 
the surface and for two different values of the bending rigidity k. The errors quoted 
in the table are due to statistical fluctuations and are estimated using the binning 
method. In Figure Q we plot a for k = 1.1. We see marked changes in a as the 
lattice size increases and as the bending rigidity increases. These discrepancies are 
a measure of the uncertainty in our determination of the Poisson ratio. We remark 
that theoretical arguments |10|ll4j indicate that the behavior in the whole flat phase 
is governed by an infrared-stable fixed point at k = 00 (the fluctuations stiffen the 
surface at long wavelength). Hence we expect to find the correct asymptotic behavior 
as L — > 00 anywhere above the crumpling transition. Larger values of the bending 
rigidity influence the results since the auto-correlation times are longer — it is more 
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Figure 1: The Poisson ratio for re = 1.1 and d = 3L/4. 



time consuming to gather the equivalent number of independent configuration. For 
both values of k the results are consistent with an approach to a value of a ~ —0.32, 
which agrees with the theoretical prediction of — 1/3 (LR) . The infinite volume value 
is extracted from the fit to <j(L) = <r(oo) + a{n)/L b . The statistical errors on the 
data do not allow for a reliable estimate of the exponent b. 

Zhang, Davis and Kroll have obtained a different result a = —0.15(1) for the 
Poisson ratio of a tethered membrane in a molecular dynamics simulation with 
periodic boundary conditions 18 . There are several possible explanations for the 
discrepancy. We stress that the present simulation goes to much larger lattice sizes: 
in view of the finite size effects demonstrated in Table [HJ we believe size to be 
relevant. Furthermore, it is hard to compare the value of the bare parameters used 
in the two simulations because both the numerical techniques and the models studied 
differ. 

For completeness we mention that Boal, Seifert and Shillcock have investigated 
numerically, and with mean field techniques, the Poisson ratio of 2-dimensional 
networks under tension |T2] . These models differ from the one we study in that they 
are strictly planar. It was found, nonetheless, that entropic effects drive the Poisson 
ratio negative. 

A comment now on our treatment of boundary effects. As noted by Abraham 
|19j . large edge fluctuations might influence the results, even though boundary effects 
should vanish in the infinite volume limit. Zhang, Davis and Kroll simulate a surface 
in a closed cell [IS], imposing periodic boundary conditions in the x-y plane only, 
and, in some cases, dynamically modifying the cell size. In our simulations the 
surface's boundary is free to fluctuate, and we need to carefully analyze the data 
in order to quantify the effect of boundary fluctuations. To this end, we restrict 
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Figure 2: The Poisson ratio for the 14 subsets at k = 1.1 and L = 46. The circles 
correspond to definition Eq. (JHJ), while the squares to definition Eq. (fTU|) . 



the definition of the correlation functions in Eq. (JHJ) to a hexagonal subset of the 
mesh (the method is described in great detail in ref. [2]). The subset excludes the 
portions of the surface close to the boundary. For each size simulated (L = 32, 
46, 64, and 128) we construct 14 hexagonal regions of increasing diameter centered 
with respect to the bulk lattice. We then measure the correlation functions (and 
hence a) restricted to these subsets. We found that, within our statistical error, 
the Poisson ratio for a subset of diameter d is in fair agreement with the one of an 
independent surface of size L = d. Thus we conclude that a does not suffer much 
from boundary effects: presumably these effects cancel out in taking the ratio of 
the correlation functions. The values quoted in Table ITT1 are measured on subsets of 
diameter d = 3L/4. The justification for this choice of subset size is given in [Oj. We 
show there that the scaling exponents are insensitive to d for subsets with d ~ 3L/4. 
Our choice is a compromise between excluding the boundary and including enough 
bulk surface for self-averaging. In Figure [21 we show the Poisson ratio for L = 46 
measured on the 14 subsets. 

Details of the statistics gathered for the various lattice sizes and different bending 
rigidities are reported in |JjJ. For the largest lattices we ran 74 x 10 6 Monte Carlo 
sweeps equivalent to 150 independent configurations. One sweep corresponds to a 
Metropolis update of each node of the lattice. 
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kindly provided computational facilities. The research of GT was sponsored by 
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of MB and MF was supported by the Department of Energy U.S.A. under contract 
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